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Abstract-  Characteristic  modes  developed  by 
Garbacz,  Harrington  and  Mautz  have  long  been 
used  in  the  analysis  of  radiation  and  scattering 
from  conducting  bodies  and  apertures.  For  their 
computation,  it  is  necessary  to  solve  an  eigen- 
system  of  the  form  Ax  — XBx.  If  the  matrices 
(A,B)  are  Hermitian  and  B is  positive  definite, 
the  generalized  eigenvalue  problem  can  be  accu- 
rately solved  using  the  simultaneous  diagonal- 
ization technique  (SDT).  Because  of  numerical 
approximations  and  rounding  sometimes  it  may 
happen  that  the  matrices  properties  deteriorate 
and  the  SDT  procedure  becomes  inapplicable. 
In  this  work  a new  technique,  developed  recently 
by  Higham  and  Cheng  is  proposed  as  a method 
to  solve  these  deteriorate  cases.  It  is  applied  to 
the  computation  of  the  characteristic  modes  for 
some  scattering  problems.  Results  are  analyzed 
and  discussed. 


1.  Introduction 

The  generalized  eigenvalue  problem: 

Ax  = ABx  (1) 

is  often  encountered  in  computational  electro- 
magnetics [1].  The  solution  of  (1)  depends  on 
properties  of  the  matrices  A and  B resulting 


from  the  discretization  of  the  field  equations.  If 
these  matrices  are  not  singular,  then  (1)  can  be 
reduced  to  the  solution  of  the  standard  eigen- 
problem: 

Hx  — Xx  with  H = B_1  A (2) 

However,  if  A and  B are  Hermitian  and  B is 
positive  definite,  the  problem  (1)  can  be  solved 
with  a higher  degree  of  accuracy  applying  the  si- 
multaneous  diagonalization  technique  (SDT)  to 
the  matrices  (A,  B)  than  the  direct  calculation 
of  the  inverse  matrix  B"1  [2].  Although  the 
underlined  properties  are  sometimes  held  of  the 
matrices  issued  in  electromagnetic  problems,  un- 
fortunately the  discretization  performed  by  nu- 
merical algorithms  on  the  electromagnetic  field 
equations  often  causes  the  loss  of  these  proper- 
ties. As  a consequence,  the  SDT  becomes  in- 
applicable. Recently,  Higham  and  Cheng  have 
faced  the  above  problem  from  a theoretical  point 
of  view  and  they  have  proposed  a technique  that 
allows  to  use  the  SDT  even  in  presence  of  dis- 
cretization errors.  This  method  is  based  on  the 
concept  of  nearest  definite  pair  [3].  In  this  work, 
we  have  applied  Higham-Cheng  results  to  evalu- 
ation of  the  characteristic  modes  for  the  scatter- 
ing from  conducting  bodies.  The  paper  is  orga- 
nized as  follows:  section  II  presents  the  basics  on 
the  theory  of  characteristic  modes  for  conduct- 
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ing  bodies.  Section  III  and  IV  illustrate  the  SDT 
procedure  for  the  resolution  of  the  positive  defi- 
nite Hermitian  generalized  eigenproblem  and  its 
extension.  In  section  V we  show  the  numerical 
results.  Finally,  in  section  VI,  the  conclusions. 


2.  Theory  of  Characteristic  Modes  for 
Conducting  Bodies 

The  characteristic  modes  have  been  intensively 
used  in  the  analysis  of  radiation  and  scatter- 
ing from  conducting  bodies  and  apertures  [4]- 
[13].  They  are  numerical  entire  basis  functions 
that  in  principle  can  be  computed  for  any  ob- 
ject shape.  Since  these  eigenfunctions  include 
the  behavior  of  unknown  current  density  flow- 
ing on  the  metallic  bodies  and  apertures,  only 
a small  number  of  them  are  required  for  a good 
reconstruction  of  it.  So,  if  available,  character- 
istic modes  would  lead  to  a scattering  or  radia- 
tion problem  treatable  with  a classical  Method 
of  Moment  even  in  presence  of  a large  number  of 
objects  [10]-[13].  In  the  case  of  scattering  or  ra- 
diation from  conducting  bodies,  these  modes  are 
basically  solutions  of  a generalized  eigenproblem 
involving  the  impedance  operator  Z , which  re- 
lates the  surface  current  J on  a conducting  body 
to  the  tangential  component  of  the  incident  elec- 
tric field  Eqllcn  on  it,  i.e.  [5]: 

ZJ  = E& ' (3) 

One  method  to  solve  a functional  equation  of  the 
form  (3)  is  to  obtain  a modal  solution  in  terms 
of  eigenfunctions  of  Z . For  this  purpose,  we 
consider  the  operational  eigenproblem: 

ZJk  — Vk MJk  (4) 

in  which  M is  a suitable  weight  operator  and 
Vk  = 1 + jA*..  The  operator  Z can  be  uniquely 
represented  in  the  form: 


Note  that  71  and  X are  real  selfadjoint  operators 
and  furthermore  71  is  positive  definite  [5].  The 
choice  M — 71  reduces  the  complex  operatorial 
eigenproblem  (4)  to  a real  operatorial  eigenprob- 
lem: 

XJk  = (8) 

The  eigenfunctions  Jk  have  been  named  char- 
acteristic modes  [5].  They  can  be  numerically 
evaluated  through  the  reduction  of  the  opera- 
tor equation  (8)  to  a matrix  equation  using  the 
Method  of  Moment  [5], [10].  For  this  aim  the  k- 
th  modal  current  density  on  the  metallic  body 
is  expanded  in  a set  of  N suitable  subsectional 
basis  functions  Bn : 

JV 

J*  = ^JnBn  (9) 

n~  1 

where  the  In  are  unknown  complex  constants. 
Substitution  of  (9)  into  (8)  and  application  of 
the  Galerkin  method  leads  to  a matrix  equation 
of  the  form: 

XJfc  = XkRIk  (10) 

where  Ik  is  the  column  vector  of  the  unknown 
coefficients  In.  Equation  (10)  is  a generalized 
eigensystem  of  the  form  (1)  in  which  the  eigen- 
values A k and  the  eigenvectors  Ik  approximate 
the  eigenvalues  and  the  eigenfunctions  of  the 
equation  (8).  It  follows  by  X and  7 Z properties 
that  the  matrices  X and  R are  expected  Her- 
mitian and  furthermore  R is  expected  positive 
definite.  Consequently,  in  principle  it  is  possible 
to  use  the  SDT  for  solving  the  eigensystem  (1). 
But  because  of  numerical  approximation,  the  R 
matrix  often  becomes  indefinite  and  the  SDT  al- 
gorithm becomes  inapplicable.  In  the  following, 
it  will  be  shown  that  the  use  of  the  nearest  def- 
inite pair  technique  permits  to  employ  the  SDT 
giving  a significant  improvement  in  the  evalua- 
tion of  the  characteristic  modes  over  the  direct 
inversion  technique. 


Z = 7l  + jX  (5) 

where  the  operators  71  and  X are  defined  as: 


3.  Resolution  of  the  Positive  Definite 
Generalized  Eigenproblem  with  the 
Simultaneous  Diagonalization  Technique 


7Z  — 
X = 


z + z * 
2 

2 j 


(6)  The  generalized  eigenproblem  where  the  matri- 
ces A and  B are  Hermitian  and  in  which  one 
is  positive  definite  plays  an  important  role  in 

(7) 

matrix  theory.  If  we  assume,  without  loss  of 
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generality,  that  B is  positive  definite,  the  prob- 
lem (1)  can  be  accurately  solved  by  means  of 
a method  named  simultaneous  diagonalization 
technique  [2]  as  follows:  Let  Ui  be  an  unitary 
matrix  whose  columns  are  an  orthonormal  set 
of  eigenvectors  for  B.  Premultiplying  and  post- 
multiplying  (1)  by  UJ  and  Ui,  where  the  star 
denotes  the  complex  conjugate  transpose  opera- 
tion, and  taking  into  account  that: 

UJUi  = UiUJ  - diag{\)  (11) 

where  diag(  1)  is  the  identity  matrix,  we  obtain: 

Ay  = AB  y (12) 

in  which  the  matrices  A,  B and  the  vector  y are 
of  the  form: 

A^UJAUi 

V = U lx  (13) 

B - UJBUi  = diagfa) 

and  in  which  diag{fi)  is  the  diagonal  matrix  formed 
by  eigenvalues  of  the  B matrix.  Next,  we  intro- 
duce the  nonsingular  transformation  H: 

H — diag(fjT^)  (14) 

Using  H in  (12),  we  obtain: 

A z = XBz  (15) 

where  A,  z and  B are: 

A = H*AH 

i = H ~ly  (16) 

B = H*BH  = diag(  1) 

Finally,  if  we  construct  an  unitary  matrix  U2 
whose  columns  are  an  orthonormal  set  of  eigen- 
vectors for  A,  the  matrix  transformation  T = 
U1HU2  permits  simultaneously  to  reduce  the 
matrices  A,  B to  diagonal  form: 

T*AT  = diag(X)  \ 

T*BT  = diag{\)  j ( ' 

resolving  the  eigenproblem  (1).  An  efficient  im- 
plementation of  the  above  procedure  that  uti- 
lizes both  the  Cholesky  factorization  and  the 
symmetric  QR  algorithm  is  described  in  [14]. 

4.  Extension  of  the  Simultaneous 
Diagonalization  Technique 

The  application  of  the  procedure  discussed  in 
the  previous  section  is  not  limited  to  the  treat- 
ment of  Hermitian  pairs  (A,B)  with  B positive 


definite,  but  sometimes  it  can  be  extended  to  the 
Hermitian  pairs  in  which  the  B matrix  is  indef- 
inite. For  this  aim,  the  concept  of  definiteness 
has  been  restated  referring  directly  to  the  ma- 
trix pair  (A,  B)  rather  than  to  a single  matrix 
B.  More  in  details,  pair  (A,  B)  having  Crawford 
number  7 [15],  denoted  as: 

7(A,B)—  min  \/ (z* Az)2  + (z*  B z)2  (18) 
IM|2=i 


strictly  positive  has  been  named  as  definite  and, 
when  7 = 0 the  pair  has  been  named  as  indefi- 
nite. Also,  it  is  possible  evaluate  the  definiteness 
of  a Hermitian  pair  (A,B)  drawing  its  Field  of 
Values  F(A  4-  jB)  in  the  complex  plane.  The 
field  of  values  is  defined  as  the  set  of  all  the  val- 
ues assumed  by  the  Rayleigh  quotients  of  the 
pair  [3]: 

F(A+  jB)  = z±  0 (19) 

z*z 

The  pair  will  be  definite  if  and  only  if  zero  does 
not  lie  in  F(A  -fjB),  otherwise  it  will  be  indef- 
inite. A solution  of  (1)  in  which  B is  indefinite, 
but  the  Crawford  number  7 of  the  pair  (A,  B)  is 
positive,  has  been  given  in  [15]  by  means  of  the 
Stewart’s  theorem  [16].  This  theorem  ensured 
the  existence  of  an  angle  0 € [0, 2zr)  that  per- 
mits to  define,  starting  from  a given  pair  (A,  B) 
with  7 > 0,  a new  pair  (A#,  Bp)  in  which  Bp  is 
positive  definite,  as  follows: 


Ap  = Acosf?  + Bs\n6  1 . . 

Bp  = - Asin0  + Bcosfl  J ^ 

For  this  new  pair  is  now  applicable  the  SDT. 
The  eigenvalues  A of  the  old  pair  (A,B)  and 
the  eigenvalues  Ap  of  the  new  pair  (Ap,Bp)  are 
simply  related  by: 


Apcos#  — sin0 
Apsin#  ~h  cos# 


(21) 


However,  it  may  happen  that  the  discretization 
errors  may  affect  the  matrices  in  a way  that  a 
pair  (A,B)  theoretically  expected  to  be  defi- 
nite becomes  indefinite  in  practice.  As  a con- 
sequence the  SDT  cannot  be  applied.  For  this 
last  case  it  can  be  useful  evaluate  the  nearest 
definite  pair  (A  -f  AA,  B + AB)  and  attempt 
to  resolve  the  problem  (1)  for  this  one.  Else, 
assigned  the  indefinite  Hermitian  pair  (A,B), 
we  must  find  the  nearest  Hermitian  definite  pair 
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(A  + AA,  B + AB),  having  a specified  Crawford 
number  7 = 8 > 0 according  to  the  distance 

[3], [17]: 

d$(A,B)  — min{\\[AA  AB]||2  : 

7(A  + AA,B  + AB)  = <5}  (22) 

and  the  angle  0 by  means  of  which  the  nearest 
pair  can  be  simultaneously  diagonalized.  This 
problem  has  been  recently  solved  by  Higham  and 
Cheng  by  means  of  the  following  theorem  [3] : 

Higham-Cheng’s  Theorem:  Let  A,B  € Cnxn 
be  Hermitian  and  let  C = A + jB  and  A<f>  = 
Aco${<p)  + Bsin(<j>).  Let  min(0<^<27r)  ^ max  (A  *) 
be  attained  at  the  angle  0 and  let  A#  have  the 
spectral  decomposition: 

A$  = Qdiag(fik)Q * with  fxn  < ...  < fii  (23) 

If  0 € F( C)  then: 

d*(A,B)  = <5  + Ml  (24) 

If  0 £ F( C)  then: 

ds  ( A , B)  = max  ( 6 + /xi , 0)  (25) 

In  both  cases,  two  sets  of  optimal  perturbations 
in  (22)  are: 

AAi  = Q diag(min(—6  -f  /x*,  O))Q*cos0  1 
ABi  = Qdiag(min{—6  + /Xfc,0))Q*sin#  J 

(26) 

and 

AA2  = -ds(A,B)diag(l)cos0  \ , . 

AB2  = —ds(A,B)diag(l)smO  J K4() 

A detailed  demonstration  of  this  theorem  can  be 
found  in  [3].  The  most  important  consequence 
of  the  above  results  is  that  the  SDT  can  be  gen- 
eralized to  the  Hermitian  pairs  (A,  B)  having  B 
matrix  indefinite. 

5.  Numerical  Results 

In  order  to  check  the  validity  and  accuracy  of  the 
foregoing  procedure  in  the  computation  of  char- 
acteristic modes  for  conducting  bodies,  numeri- 
cal results  for  some  cases  are  carried  out.  Firstly, 
characteristic  modes  have  been  computed  for  the 
TMZ  scattering  from  a perfectly  conducting  cir- 
cular cylinder  and  from  an  elliptic  one,  since  for 


Figure  1:  Geometry  for  scattering  from  circular, 
elliptic  and  square  metallic  cylinders. 


Figure  2:  Plot  of  the  smaller  normalized  eigenvalues 
for  the  TMZ  scattering  from  a circular  cylinder  of 
radius  a = ^ (Problem  Size  128  x 128), 

these  cases  characteristic  modes  are  known  in 
analytical  form  [12].  In  Fig.(l)  are  represented 
the  geometries  for  both  problems.  In  Fig. (2)  is 
reported  the  plot  of  the  smaller  eigenvalues  of 
the  R matrix  obtained  discretizing  the  EFIE  for 
the  scattering  due  to  a plane  wave  impinging 
on  a metallic  circular  cylinder  with  radius 
Pulse  basis  functions  and  delta  function  test- 
ing functions  sire  used  to  create  the  matrix  R 
[1].  Since  only  the  properties  of  the  matrices 
above  are  considered  here  for  the  computation 
of  characteristic  modes,  it  is  not  necessary  to 
specify  the  direction  of  the  incident  field.  The 
depicted  eigenvalues  are  scaled  to  show  the  rel- 
ative value  only.  This  plot  clearly  shows  that 
the  R matrix  is  indefinite,  even  if  it  is  expected 
definite  because  the  properties  of  the  operator 
H.  In  Fig.  (3)  is  plotted  the  field  of  values  for 
the  pair  (X,R).  Because  the  origin  of  the  com- 
plex plane  is  not  contained  in  F(X  + jR)  this 
pair  is  definite.  In  Figs. (4)  and  (5)  are  reported 
the  plot  of  the  smaller  eigenvalues  of  the  R ma- 
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Figure  3:  Field  of  Values  for  pair  related  to  compu- 
tation of  characteristic  modes  for  circular  cylinder. 

trix  and  the  field  of  values  F(X  + jH)  for  the 
pair  obtained  considering  the  scattering  from  a 
metallic  elliptic  cylinder  with  semi  major  axis 
a = A and  semi  minor  axis  b = ^ by  an  inci- 
dent plane  wave.  For  this  scattering  problem  we 
have  obtained  both  an  indefinite  matrix  B and 
an  indefinite  pair  (X,R)  as  clearly  depicted  in 
these  figures.  In  Table  (1)  are  shown  the  an- 
alytical eigenvalues  versus  the  numerical  ones 
for  the  previously  analyzed  scattering  problems. 
Numerical  eigenvalues  have  been  obtained  using 
both  the  Higham-Cheng  procedure  and  the  di- 
rect inversion  technique.  Excellent  agreement  is 
shown  comparing  the  exact  eigenvalues  and  the 
numerical  ones  computed  by  means  of  the  former 
procedure  while  the  eigenvalues  obtained  using 
the  direct  inversion  technique  are  very  inaccu- 
rate. As  last  case,  consider  the  evaluation  of 
the  characteristic  modes  for  the  scattering  from 
a metallic  square  cylinder  of  side  a = A.  In 
Figs.(l)  and  (6)  are  reported  the  geometry  for 
the  problem  and  its  field  of  values,  respectively. 
As  it  is  shown  in  Fig. (6),  the  matrix  pair  (X,R) 
for  this  case  is  definite  because  the  field  of  val- 
ues does  not  contain  the  origin  of  the  complex 
plane.  In  Fig.  (7)  is  shown  the  obtained  current 
distribution  on  the  cylinder  by  a plane  wave  im- 
pinging at  an  angle  ==  0.  The  continuous  line 
indicates  the  current  computed  by  means  of  the 
standard  MoM  while  crosses  indicate  the  solu- 
tion obtained  using  characteristic  modes  com- 
puted using  the  Higham-Cheng  procedure.  A 
very  good  agreement  is  clearly  observed. 
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Figure  4:  Plot  of  the  smaller  normalized  eigenvalues 
for  the  T Mz  scattering  from  an  elliptic  cylinder  with 
semi  major  axis  a = A and  semi  minor  axis  6 = £ 
(Problem  Size  180  X 180). 


Figure  5:  Field  of  Values  for  pair  related  to  compu- 
tation of  characteristic  modes  for  elliptic  cylinder. 

6.  Conclusions 

A technique  developed  by  Higham  and  Cheng 
for  the  treatment  of  the  Hermitian  generalized 
eigenproblem  based  on  the  concept  of  nearest 
definite  pair  has  been  proposed  for  the  compu- 
tation of  the  characteristic  modes.  Numerical  re- 
sults for  some  scattering  problems  are  presented. 
In  all  the  cases,  the  outlined  procedure,  even 
if  acting  on  very  ill  conditioned  matrices,  pro- 
vides very  accurate  numerical  results.  Finally, 
we  point  out  that  the  application  of  the  method 
is  not  limited  to  the  presented  cases  being  fruit- 
fully applicable  to  a wide  number  of  electromag- 
netic problems  which  posses  the  discussed  prop- 
erties [18]. 
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Figure  6:  Field  of  Values  for  pair  related  to  com- 
putation of  characteristic  modes  for  square  cylinder 
(Problem  Size  165  x 165). 


Figure  7:  Plot  of  the  current  density  on  the  metallic 
square  cylinder.  The  continuous  line  indicate  solu- 
tion computed  by  means  of  the  MoM.  Crosses  indi- 
cate solution  obtained  using  characteristic  modes. 
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-2.58841 

Table  1:  Analytical  vs.  numerical  eigenvalues. 
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